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Abstract 


Traditionally, stellar radiation has been the only heat source considered capable of determining global climate on 
long timescales. Here, we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high- 
enough levels to induce a runaway greenhouse for a long-enough duration for all the hydrogen to escape. 
Without hydrogen, the planet no longer has water and cannot support life. We call these planets “Tidal Venuses” 
and the phenomenon a “tidal greenhouse.” Tidal effects also circularize the orbit, which decreases tidal heating. 
Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their 
water, but eventually settle into nearly circular orbits (i.e., with negligible tidal heating) in the habitable zone 
(HZ). However, these planets are not habitable, as past tidal heating desiccated them, and hence should not be 
ranked highly for detailed follow-up observations aimed at detecting biosignatures. We simulated the evolution 
of hypothetical planetary systems in a quasi-continuous parameter distribution and found that we could con- 
strain the history of the system by statistical arguments. Planets orbiting stars with masses <0.3 Msun may be in 
danger of desiccation via tidal heating. We have applied these concepts to Gl 667C c, a ~4.5 Megarin planet 
orbiting a 0.3 Msun star at 0.12 AU. We found that it probably did not lose its water via tidal heating, as orbital 
stability is unlikely for the high eccentricities required for the tidal greenhouse. As the inner edge of the HZ is 
defined by the onset of a runaway or moist greenhouse powered by radiation, our results represent a funda- 
mental revision to the HZ for noncircular orbits. In the appendices we review (a) the moist and runaway 
greenhouses, (b) hydrogen escape, (c) stellar mass-radius and mass-luminosity relations, (d) terrestrial planet 
mass-radius relations, and (e) linear tidal theories. Key Words: Extrasolar terrestrial planets—Habitability— 
Habitable zone—Liquid water—Tides. Astrobiology 13, 225-250. 


1. Introduction 


A" LIFE ON EARTH requires liquid water. Not surprisingly, 
our search for life on exoplanets therefore begins with 
identifying environments that support it (Hart, 1979; Kasting 
et al., 1993; Selsis et al., 2007; von Bloh et al., 2007). Many other 
features of a planet are also important, such as the mix of gases 
in the atmosphere and the interior’s structure and energy 
budget, but current research has focused primarily on the 
stability of surface water. On Earth, liquid water persists 
primarily because the Sun’s radiation heats the surface to a 
temperature between water’s freezing and boiling points. 
Thus, the concept of a “habitable zone” (HZ) emerged, which 
is the region around a star in which insolation can maintain 
liquid water on the surface, assuming an Earth-like planet. 


Stars have a wide range of luminosities, most of which are 
considerably lower than that of the Sun, yielding HZs that 
are closer in. At the extreme, the luminosity becomes so low 
that a planet orbiting at the stellar surface would be too cold 
to support liquid water on its surface. Such a star has no HZ. 
The recently discovered Y dwarf (which is not a star) WISEP 
J1828+2650 (Cushing et al., 2011), with an effective temper- 
ature below 300 K, is an example of such a primary. For 
warmer stars, HZs may still be close enough in that non- 
radiative processes may impact habitability, such as stellar 
flaring (e.g., Khodachenko et al., 2007; Lammer et al., 2007, 
2010; Tian, 2009; Segura et al., 2010), decreased initial volatile 
inventory (Lissauer, 2007; Raymond et al., 2007), or tidal ef- 
fects (e.¢., Kasting et al., 1993; Joshi et al., 1997; Correia et al., 
2008; Jackson et al., 2008a; Barnes et al., 2009a; Heller et al., 
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2011). As the HZ of our Sun is too distant for these phe- 
nomena to affect Earth, we can currently only explore their 
role theoretically; consequently, many scientists consider 
close-in planets less favorable candidates for habitability. 

Nevertheless, the last few years have seen renewed in- 
terest in the potential habitability of planets in orbit about 
low-luminosity objects (Scalo et al., 2007; Tarter et al., 2007; 
Lunine et al., 2008; Monteiro, 2010; Agol, 2011; Bolmont et al., 
2011). This shift occurred because Earth-sized planets are 
easiest to detect around low-luminosity hosts, due to the 
larger mass and radius of the planet relative to the star. 
Furthermore, these objects are the most abundant in the solar 
neighborhood. Planets with masses between 1 and 10 Meartn, 
“super-Earths,” have indeed been detected in the last few 
years around low-mass stars, such as G1 581 d and Gl 667C c 
by radial velocity (Udry et al., 2007; Mayor et al., 2009; Vogt 
et al., 2010; Bonfils et al., 2011; Forveille et al., 2011; Anglada- 
Escudé et al., 2012; Delfosse et al., 2012) and GJ 1214 b by 
transit (Charbonneau et al., 2009; Sada et al., 2010; Carter 
et al., 2011; Kundurthy et al., 2011). Several observational 
campaigns designed specifically to detect planets around 
low-mass stars are now underway (Nutzman and Char- 
bonneau, 2008; Boss et al., 2009; Zechmeister et al., 2009; Bean 
et al., 2010; Rodler et al., 2011). 

The inner edge of the HZ is especially prone to non- 
radiative phenomena. Traditionally, the inner boundary is 
defined by the stellar distance at which insolation is strong 
enough to remove all water in the atmosphere and on the 
surface (Kasting, 1988). Two different dehydrating scenarios 
(the “moist” and “runaway” greenhouse) are discussed in 
detail in Section 2 and Appendix A. Both require water vapor 
to penetrate a stratosphere, be dissociated (photolyzed) by 
high-energy radiation, and culminate in the escape of hy- 
drogen. We call any process that can ultimately lead to total 
water loss a “desiccating greenhouse.” Without sufficient 
hydrogen, water cannot form, and the planet will remain 
uninhabitable forever, unless a major, unlikely event occurs, 
for example, an impact by a water-rich body that simulta- 
neously delivers water and changes the atmosphere in a way 
that halts the desiccating greenhouse. This definition of the 
inner edge is conservative because no known process can 
maintain habitability against a desiccating greenhouse. Other 
processes could be equally deleterious for life, but clearly 
total desiccation will terminate habitability. 

We will discuss four types of terrestrial planets in this 
study, classified by their water content. “Wet” planets are 
terrestrial exoplanets that have a water content similar to that 
of Earth. “Dry” exoplanets have far less, ~3cm deep if 
condensed and spread globally, but are habitable (see Abe 
et al., 2011). “Desiccated” planets have lost all their water 
through a desiccating greenhouse. Finally, “water worlds,” 
(Léger et al., 2004; Raymond et al., 2004) are planets with a 
much larger inventory of water than modern Earth, for ex- 
ample, a warm Europa. This investigation explores the 
transition of wet, dry, and water worlds into the desiccated 
State. 

Dole (1964) was the first to point out that terrestrial 
planets in the HZ of low-luminosity stars can have their spin 
altered by tidal interaction. In particular, the danger of 
synchronous rotation, that is, one hemisphere always facing 
the star, was emphasized. Kasting et al. (1993) quantified this 
concept and found that planets orbiting within the HZ of 
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stars less than two-thirds the mass of the Sun were in danger 
of synchronization. Although their analysis was limited to 
Earth-like planets on circular orbits, a general belief devel- 
oped that those planets could not be habitable, as one half of 
the planet would freeze while the other would roast. 

Synchronization is certainly an important consideration 
when assessing habitability, and many investigations have 
explored its role, but with mixed results. Atmospheric 
modeling initially suggested that circulation will transport 
energy to the unlit side, ameliorating the extreme tempera- 
ture difference (Joshi et al., 1997). Some subsequent modeling 
has confirmed that synchronous rotators are likely to have 
super-rotating atmospheres (Joshi, 2003; Edson et al., 2011; 
Heng and Vogt, 2011; Showman and Polvani, 2011), while 
others have discounted it (Wordsworth ef al., 2011), and still 
others have suggested that such a state might be beneficial at 
the outer edge of the HZ (Pierrehumbert, 2011). Taken to- 
gether, these investigations suggest that synchronized plan- 
ets should not be dismissed uniformly as uninhabitable. 
Hence, synchronization is not as stringent a constraint as 
desiccation and therefore is not an HZ boundary. 

For many years, confusion also existed regarding the term 
“tidal locking.” Many investigators assumed it was synon- 
ymous with “synchronous rotation.” If the orbits are non- 
circular, as for many exoplanets (Butler et al., 2006), then 
tidally evolved planets may reach an equilibrium state where 
they rotate faster than synchronous with an “equilibrium” or 
“pseudo-synchronous” period. This aspect of tidal theory 
has been known for decades (e.g., Goldreich, 1966; Greenberg 
and Weidenschilling, 1984) but has only recently been 
pointed out for the case of exoplanets (Barnes et al., 2008; 
Correia et al., 2008; Ferraz-Mello et al., 2008). Therefore, some 
exoplanets, such as Gl 581 d with an eccentricity of 0.38 
(Mayor et al., 2009), may be “tidally locked” but rotate about 
twice per orbit (Barnes et al., 2008; Heller et al., 2011). For 
more on this point, see the work of Murray and Dermott 
(1999), Chapter 5.2; and for a counterargument, see the work 
of Makarov and Efroimsky (2012). In the present paper, we 
use “tidally locked” to mean a planet rotating at the equi- 
librium period as determined by its eccentricity and obliq- 
uity; see Eqs. E7, E15-E16. In summary, a synchronously 
rotating planet is tidally locked (yet it could still have non- 
zero eccentricity and obliquity), but a tidally locked planet is 
not necessarily rotating synchronously. 

Even if an orbit is currently circular, tides may not have 
driven the rotation to synchronous. If the orbit began with 
large eccentricity, tides will tend to damp it to zero, and we 
may expect it to be rotating synchronously. However, the 
planet could pass through one or more “spin-orbit reso- 
nances,” where the planet’s rotational frequency is com- 
mensurate with its orbital frequency (see, e.g., Rodriguez 
et al., 2012). For example, Mercury rotates three times for 
every two times it orbits the Sun, a 3:2 spin-orbit resonance. 
Spin-orbit resonances require a radially asymmetric mass 
distribution, which is likely for tidally deformed exoplanets 
but cannot be measured for the foreseeable future. A planet 
caught in a spin-orbit resonance may remain in that state 
even if circularized, as a resonance is a strong dynamical 
process. For any particular exoplanet, capture and retention 
into a spin-orbit resonance will be very difficult to constrain 
observationally, so all reasonable options should be consid- 
ered. Therefore, synchronous rotation is unlikely for planets 
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with large eccentricities (Mercury’s is 0.2) and not even 
guaranteed for a circular orbit. For more on spin-orbit reso- 
nances, the reader is referred to the work of Murray and 
Dermott (1999), Chapter 5.4. 

As tidal locking of the planetary rotation is not an absolute 
constraint on habitability, we turn to tidal heating as the 
other tidal phenomenon most likely to affect planetary hab- 
itability. As a planet moves from periastron, its closest ap- 
proach to the star, to apoastron, the farthest point, and back 
again, the gravitational force changes, being inversely pro- 
portional to distance squared. This difference creates an os- 
cillating strain on the planet that causes its shape to vary 
periodically. The rigidity of the planet resists the changes in 
shape, and friction generates heat. This energy production is 
called tidal heating. 

Tidal heating is responsible for the volcanism on Io (Strom 
et al., 1979; Laver et al., 2007), which was predicted by Peale 
et al. (1979) with the use of tidal theory. Io is a small body 
orbiting Jupiter with an eccentricity of 0.0041 that is main- 
tained by the gravitational perturbations of its fellow Gali- 
lean moons and shows global volcanism that resurfaces the 
planet on a timescale of 100 to 10° years (Johnson et al., 1979; 
Blaney et al., 1995; McEwen et al., 2004). The masses of Jupiter 
and Io are orders of magnitude smaller than a star and ter- 
restrial exoplanet; thus the latter have a much larger reser- 
voir of orbital and rotational energy available for tidal 
heating. Moreover, some exoplanets have been found with 
orbital eccentricities larger than 0.9 (Naef et al., 2001; Jones 
et al., 2006; Tamuz et al., 2008). Thus, the tidal heating of 
terrestrial exoplanets may be much more effective than that 
which occurs on Io (Jackson et al., 2008a, 2008c; Barnes et al., 
2009a, 2010; Heller et al., 2011). This expectation led to the 
proposition that terrestrial exoplanets with surface heat 
fluxes as large or larger than Io’s should be classified as 
“Super-los” rather than “Super-Earths” (Barnes et al., 2009b). 
Numerous Super-Io candidates exist, such as CoRoT-7 b 
(Léger et al., 2009; Barnes et al., 2010), Gl 581 e (Mayor et al., 
2009), 55 Cne e (McArthur et al., 2004; Dawson and Fab- 
rycky, 2010; Winn et al., 2011), and Kepler-10 b (Batalha et al., 
2011), but none is in the HZ. 

Jackson et al. (2008a) and Barnes et al. (2009a) considered 
Io’s heat flux, ~2 W m~? (Veeder et al., 1994; Spencer et al., 
2000; McEwen et al., 2004), to be an upper limit for habit- 
ability, arguing that Io-like surfaces are dangerous for hab- 
itability. However, 2 W m * may not be sufficient to sterilize 
a planet and so should not be considered a hard limit to 
habitability. For example, the heat flow in inhabited hy- 
drothermal vent systems on Earth, such as the Endeavour 
segment of the Juan de Fuca Ridge (Holden et al., 1998), is 
~30 W m ° (Fontaine et al., 2011). Thus, a water world with 
tens of watts per square meter of energy output could 
support life. While Io-like volcanism is clearly an issue for 
habitability, it may not always lead to sterilization. How- 
ever, if the tidal heating can maintain a desiccating green- 
house long enough for all the water to be lost, then the 
planet becomes uninhabitable and is highly unlikely to ever 
regain habitability. 

Calculations of the inner edge of the HZ have traditionally 
assumed the primary energy source at the surface is stellar 
radiation, as is the case for Earth. Following Barnes et al. 
(2009a), we call that type of HZ an “insolation habitable 
zone” (IHZ), as starlight is the only energy source consid- 
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ered. In this investigation, we identify the amount of tidal 
heating that triggers a desiccating greenhouse, as well as the 
combinations of physical and orbital parameters for which 
tidal heating could yield it. We find such a state is predicted 
by current models and dub such a world a “Tidal Venus” 
and the phenomenon that produces it a “tidal greenhouse.” 
The tidal greenhouse would probably have the same effect 
on habitability as one caused by irradiation and hence 
should be considered as hard a limit to habitability as a 
traditional, insolation-driven desiccating greenhouse. 

In this study, we define the limits of Tidal Venuses in 
terms of stellar and planetary mass, M,, and Mp, respectively, 
the orbital semimajor axis a, orbital eccentricity e, planetary 
radius R,, planetary obliquity wp, and planetary spin fre- 
quency W,. Some of these quantities are easily observed by 
current technology; others require special geometries and the 
next generation of space telescopes. Therefore, for the next 
decade, not all newly found planets in an IHZ will have well- 
constrained tidal heating. This study provides a framework 
for identifying the range of tidal heating on terrestrial planets 
orbiting in the IHZ of low-luminosity stars. As can be seen 
below, the story is complicated and involves a large pa- 
rameter space, but is tractable. 

We first review the surface conditions that lead to a 
greenhouse state in Section 2, including an extended dis- 
cussion of the moist and runaway greenhouses (Appendix 
A); a short review of hydrogen escape driven by high-energy 
radiation (Appendix B); the relationships between mass, ra- 
dius, and luminosity for low-mass, hydrogen-burning stars, 
“M dwarfs,” or “late-type stars” in the parlance of astro- 
physics; the extent of the [HZ (Appendix C); and mass- 
radius relationships for terrestrial exoplanets (Appendix D). 
We then briefly describe tidal heating in Section 3 with de- 
tails in Appendix E. Next, we show that tidal heating alone 
can produce surface conditions capable of triggering a run- 
away greenhouse on planets orbiting M dwarfs (Section 4). 
Then we explore how past tidal heating may preclude hab- 
itability of planets found in the IHZ regardless of their cur- 
rent tidal heating (Section 5). In Section 6, we then consider 
the Gl 667C system, which contains two potentially habitable 
planets, and find that tidal heating is unlikely to have ster- 
ilized either. In Section 7, we discuss the results, and finally 
in Section 8 we draw our conclusions. 


2. The Inner Edge of the Habitable Zone 


Planets with surface water will always have water vapor 
in their atmospheres. The amount of water vapor present in 
equilibrium with a liquid surface is described by the satu- 
ration vapor pressure P.ar, Which is the pressure exerted by 
water vapor in thermodynamic equilibrium with standing 
surface water. The value of psa; depends exponentially on 
temperature: a warmer planet will have much more water 
vapor in its atmosphere. Water vapor is a greenhouse gas 
(see below), so a warmer planet will have a stronger green- 
house effect, which enhances warming, a positive (but not 
necessarily runaway) feedback. 

Atmospheric gases on Earth—and those that we expect on 
habitable Earth-like planets—are mostly transparent to op- 
tical wavelengths, so starlight is able to heat the surface. The 
surface in turn heats the air in contact with it, which rises, 
expanding and cooling adiabatically. As the air cools, water 
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vapor condenses, releasing latent heat and so slowing the 
cooling. The net rate of temperature decrease with height is 
called the “moist adiabatic lapse rate,” setting the mean 
temperature—pressure (T-p) structure of Earth’s troposphere 
(the lower region of the atmosphere affected by convection 
and bounded at the bottom and top by the planetary surface 
and the tropopause, the altitude at which temperature stops 
dropping with height). 

Atmospheric gases that absorb radiation at similar wave- 
lengths to the radiative emission of the planet (thermal IR for 
Earth and any habitable planet) are termed “greenhouse 
gases”; these cause a “greenhouse effect.” The most impor- 
tant greenhouse gases for Earth are water vapor and carbon 
dioxide. They absorb radiation emitted by the surface and 
then reemit, both downward toward the surface and upward 
into space. Energy radiated toward the surface (commonly 
called back radiation) heats the surface. Most importantly, as 
the atmosphere is cooler than the surface, the amount of 
radiation emitted to space from the atmosphere is less than 
the amount of energy emitted by the surface. Thus, the 
presence of a greenhouse atmosphere means the surface 
temperature is warmer than the equilibrium temperature of 
the planet. 

Earth’s greenhouse effect keeps the planet warm enough 
to be habitable; without it Earth’s surface would be at the 
planet’s equilibrium temperature, a barren 255 K. However, 
stronger and more water-vapor-rich greenhouse atmo- 
spheres can render the planet uninhabitable through high- 
temperature sterilization and/or desiccation (loss of the 
ocean), and thus define the inner edge of the HZ. There are 
two physically distinct situations that lead to loss of habit- 
ability, the “runaway greenhouse” and “moist greenhouse.” 

The runaway greenhouse was recently reviewed by 
Goldblatt and Watson (2012), so we provide a summary 
description only. As the planet warms, the amount of water 
vapor in the atmosphere increases such that water becomes a 
major constituent of the atmosphere and ultimately the 
dominant one. Consequently, the moist adiabatic lapse rate 
tends toward the saturation vapor pressure curve for water, 
and the T-p structure of the atmosphere becomes fixed. 
Concurrently, the atmosphere becomes optically thick in the 
thermal IR, such that only the upper troposphere can emit to 
space. As the T-p structure is fixed, the emitted radiation is 
also fixed, imposing a limit on the outgoing radiation (Fait) 
from the troposphere. Values for F,,i¢ from the literature are 
typically 285 W m~* to 310 W m ® for a 1 Megartn planet 
(Pollack, 1971; Watson et al., 1984; Abe and Matsui, 1988; 
Kasting, 1988; Pierrehumbert, 2010). The physics of this limit 
were described in more detail by Simpson (1927) and Na- 
kajima et al. (1992). [Komabayashi (1967) and Ingersoll (1969) 
described a stratospheric limit at a higher flux, 385 W in 
which is never reached in practice.] If the amount of energy 
supplied to the atmosphere by the Sun (Simpson, 1927), 
impacts (Abe and Matsui, 1988), or tidal heating (this work) 
were to exceed F.it, then the atmosphere would not be able 
to maintain radiation balance, and runaway heating of the 
surface would ensue, causing evaporation of the entire 
ocean. Radiation balance would be regained either (a) when 
the surface temperature reaches ~ 1400 K, at which point 
enough radiation is emitted in the near-IR where water va- 
por is not a good absorber or (b) if all the water vapor is lost 
from the atmosphere. 
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One likely water loss mechanism is hydrogen escape to 
space. Today, little water vapor reaches the upper atmo- 
sphere because the tropopause acts as a “cold trap”; almost 
all water has condensed lower in the atmosphere, so the 
tropopause water vapor mixing ratio (the abundance of 
water vapor in the atmosphere, expressed as the ratio of the 
mass of water vapor to the mass of dry air) is low. Water 
vapor transport above here is generally diffusive, leading to 
a constant water vapor mixing ratio in the upper atmo- 
sphere. (More precisely, in Earth’s atmosphere, water vapor 
increases with altitude in the stratosphere due to methane 
oxidation and decreases with altitude in the next higher at- 
mospheric layer, the mesosphere, due to photolysis; how- 
ever, the total hydrogen mixing ratio is conserved, and the 
hydrogen escape rate depends on this value.) Thus, water is 
not a strong source of escaping hydrogen, and the overall 
escape rate is low. However, in the runaway greenhouse the 
atmosphere becomes predominantly water, so no such con- 
straint applies, and hydrogen may escape hydrodynamically 
to space (Kasting and Pollack, 1983). The rate will depend on 
the amount of extreme ultraviolet (XUV, 1-1200 A) radiation 
absorbed in the highest part of the atmosphere, so this pro- 
cess is also limited by the available stellar XUV energy 
(Watson et al., 1981). A Venus-sized planet at the inner edge 
of the HZ could lose an ocean the size of Earth’s in ~10° 
years (Watson et al., 1981); we refer to this as the desiccation 
time, tge, (discussed further below). Deuterium would be 
retained preferentially over ordinary hydrogen during hy- 
drogen escape. Enrichment of D/H on Venus implies that 
Venus lost a substantial amount of water to space (Donahue 
et al., 1982; de Bergh, 1993), likely after experiencing a run- 
away greenhouse. 

We adopt tges= 10° years based on the work on Venus as 
well as on calculations presented in Appendix B for a ter- 
restrial planet orbiting a very low-mass star. Such stars are 
very active and emit strongly in the near UV, especially 
when they are young (e.g., West et al., 2008), but the total flux 
is still less than that of the present-day Sun (Fleming et al., 
1993), and the XUV flux (which drives hydrodynamic es- 
cape) is not well known. Thus, we expect a range of tges to 
exist, especially since the orbit evolves with time due to tides, 
but choose this value for simplicity. The model described in 
Appendix B can be applied to individual cases when obser- 
vations permit realistic modeling. We stress that there are 
many unknowns that may affect the actual value of tg... The 
masses of exoplanet oceans are unknown. Earth’s mantle 
contains a few ocean masses of water (Bell and Rossman, 
1992; Murakami et al., 2002; Marty and Yokochim, 2006), so it 
is possible that the surface of a planet could acquire a new 
ocean via outgassing after desiccation and become habitable 
again. Thus, we cannot conclude with certainty that a planet 
with a desiccated atmosphere will never be habitable in the 
future. However, we also expect rapid mantle overturn 
during strong tidal heating; hence the entire water inventory 
may actually be lost. We therefore assume that the tidal 
greenhouse is highly likely to sterilize a planet permanently. 

The moist greenhouse (Kasting, 1988) describes a warm 
atmosphere, in which the whole troposphere is assumed to 
be water vapor-saturated and is underlain by a liquid ocean. 
As the surface temperature increases, the tropopause is 
pushed higher (to lower pressure), given the reasonable 
assumption of a constant tropopause temperature. While the 
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saturation vapor pressure at the tropopause is, by this as- 
sumption, constant, the saturation mixing ratio of water at 
the tropopause (Psat/Ptrop) increases as the tropopause 
moves higher, where Prop is the atmospheric pressure at the 
tropopause. The cold trap is then no longer effective, sub- 
stantial water vapor penetrates the stratosphere, and water- 
derived hydrogen escape can be effective. The planet may 
thus gradually desiccate while in a hot, but stable, climate. 
This process could be driven by high greenhouse gas in- 
ventories rather than external heating. 

A key distinction between the runaway and moist green- 
houses is that the former happens when a known flux of 
energy is supplied to the planet, whereas the latter depends 
most strongly on surface temperature. The runaway green- 
house occurs because water vapor is a greenhouse gas and 
can therefore trap heat near the surface: hence it is only a 
function of the absorptive properties of water and the energy 
flux. The moist greenhouse occurs when the temperature at 
the tropopause is large enough that water does not con- 
dense and is therefore able to escape to the stratosphere 
where it can be photolyzed. Increasing the planet’s inven- 
tory of noncondensible greenhouse gases has a small effect 
on the runaway greenhouse limit but can drive the moist 
greenhouse since the relative amount of H2O is lower. Thus, 
the runaway greenhouse is a more conservative choice to 
demarcate the inner edge of the IHZ in the sense that it 
depends solely on water, whereas knowledge of the atmo- 
spheric composition may not be available to assess the 
likelihood of a moist greenhouse. 

In Fig. 1, we show the moist and runaway greenhouse for 
wet planets in orbit around M dwarfs. See Appendices A, C, 
and D for a discussion of IHZ limits, stellar mass-radius and 
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FIG. 1. Comparison of HZ boundaries for different plane- 
tary masses and desiccating greenhouses. The shaded regions 
represent the IHZ boundaries from Selsis et al. (2007): dark 
gray assumes no cloud coverage, medium gray 50%, and 
light gray 100%, and assuming a 1 Mga; planet. The black 
curves represent the runaway greenhouse limit from Pierre- 
humbert (2011), with a planetary albedo of 0.49. The solid 
curve is for a 30 Mgarth planet, dashed for a 0.3 Megarth planet. 
For these calculations we used the stellar mass-radius rela- 
tionship from Bayless and Orosz (2006), the mass-luminosity 
relationship from Reid and Hawley (2000), and the terrestrial 
mass-radius relationship from Sotin et al. (2007). 
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mass-luminosity relationships, and terrestrial mass-radius 
relationships. In Fig. 1, the gray regions are the limits of the 
moist greenhouse as presented in Selsis et al. (2007). At the 
inner edge, Selsis et al. (2007) found that 300 W m° triggers 
the moist greenhouse on a 1 Megartn, planet, and from left to 
right the limits assumed 100%, 50%, and 0% cloud cover (for 
a 0.25 Msgun star these limits correspond to bond albedos of 
0.75, 0.49, and 0.23, respectively). The solid curve is the 
runaway greenhouse limit (Pierrehumbert, 2010) for a 30 
Mearth planet and dashed for a 0.3 Megartn planet, both with an 
albedo of 0.49 (compare to the medium gray). As expected 
(see Appendix A), the smaller planet’s inner edge lies at 
larger semimajor axis than the 1 Mgartn moist greenhouse 
limit, which in turn lies farther out than the 30 Mg, run- 
away greenhouse limit. 

We have presented the classical description of the desic- 
cation at the inner edge of the HZ above. Various compli- 
cations are worthy of note. Water on a dry planet will get 
trapped at the poles, which makes a moist or runaway 
greenhouse harder to achieve and means that the inner edge 
of the HZ is nearer the star (Abe et al., 2011). We do not 
include that limit in Fig. 1, but we return to it in Section 4. 
Figure 1 assumes the planetary orbit is circular. However, for 
the many exoplanets with (highly) eccentric orbits (Butler 
et al., 2006), the total irradiation over an orbit determines the 
annual-averaged surface temperature (Williams and Pollard, 
2002) and pushes the IHZ boundaries out by a factor of 
(1-e7)~1/4 (Barnes et al., 2008). 


3. Tidal Heating 


Tidal heating is responsible for the volcanic activity on Io 
(Peale et al., 1979) and probably the geysers of Enceladus 
(Hansen et al., 2006; Porco et al., 2006; Hurford et al., 2007). 
Tidal theory has a long and established body of work; 
however, tidal processes remain poorly understood. The 
difficulty lies in the complexity of the energy dissipation 
processes and the very long timescales associated with tidal 
evolution. The Earth-Moon system is the most accurately 
studied, with lunar laser ranging providing a precise mea- 
surement of the recession of the Moon due to tides of 
~38mm/yr (Dickey et al., 1994), as well as direct measure- 
ments of the locations of dissipation in the ocean (Egbert and 
Ray, 2000). However, using the currently estimated tidal 
dissipation parameters to extrapolate backward in time 
predicts that the Moon was at Earth’s surface about ~2 Gyr 
ago (MacDonald, 1964), which contradicts the standard im- 
pact model for the lunar origin [note that numerous issues 
remain such as the origin of Earth’s and the Moon’s obliq- 
uities (Touma and Wisdom, 1994) and the perturbations 
from other planets (Cuk, 2007)]. The deceleration of Io’s or- 
bital velocity has also been tentatively detected and seems 
broadly consistent with tidal theory (Aksnes and Franklin, 
2001; Lainey et al., 2009). However, these data lie in the same 
regime, that of nearly circular orbits. Exoplanets have been 
found with extremely eccentric orbits, so we can only ex- 
trapolate from our Solar System cautiously. The details of 
tidal theory are complicated; hence we relegate the discus- 
sion to Appendix E. 

Two end-member models of tidal theory have been ap- 
plied to exoplanets and the bodies of our solar system: one 
assumes the energy dissipation gives a constant phase lag in 
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the periodic distortion (CPL), and the other assumes a con- 
stant time lag (CTL) (Greenberg, 2009). While simple and 
linear, such models are probably commensurate with the 
dearth of information we have about exoplanet interior 
processes. More complicated models have been constructed, 
and they reproduce the above models for certain choices of 
internal composition, structure, and energy transport (¢.g., 
Henning et al., 2009). Thus, the CPL and CTL models can 
provide important and accurate insight into the tidal evolu- 
tion of exoplanets. 

These models converge at e=0 and have been shown to be 
nearly identical for e< 0.2 (Leconte et al., 2010) and when 
using Eq. E21. However, for e>0.3, they diverge significantly. 
We urge caution when interpreting results in the upcoming 
sections, which allow e to be as large as 0.8. We include this 
range primarily for illustrative purposes and as a baseline for 
any future work that may include nonlinear effects. 

The tidal heating of a body is provided in Eqs. E10 and 
E18 for the CPL and CTL models, respectively. Averaging 
the heating rate of the entire planetary surface gives the 
surface energy flux due to tides, Fiige. In those equations, the 
strength of the tidal effects is parameterized as a “tidal 
quality factor” Q (CPL) or “tidal time lag” t (CTL), which are 
notoriously difficult to measure or estimate from first prin- 
ciples. Earth’s current values are Q=12 (Yoder, 1995) and 
t=638s (Lambeck, 1977; Neron de Surgy and Laskar, 1997), 
respectively. However, as noted above, these values predict 
too short a lifetime of the Moon. This discrepancy has led to 
the notion that Earth’s response to lunar tides has evolved 
with time, possibly due to changing size, shape, and seafloor 
topography of the oceans, affecting the ocean current’s re- 
sponse to the tidal potential. Measurements of the dry bodies 
in our solar system have found that their Q values tend to 
cluster around 100; see Appendix E. 

Recent satellite observations of Earth have revealed the 
locations of tidal dissipation in our oceans (Egbert and Ray, 
2000). Tides force water through shallow seas and straits, 
which causes energy dissipation. In the open ocean, tidal 
dissipation is probably a nonlinear process in which currents 
are disrupted by seafloor topography. The presence of the 
ocean provides more opportunity for tidal dissipation than 
on dry planets, and we assume that most dissipation on 
habitable exoplanets will also occur in their oceans. There- 
fore, for the following calculations, we assume modern 
Earth-like planets with Q=10-100 or t=64-640s. 


4. Tidal Venuses 


We computed tidal heating rates for a range of planetary, 
stellar, and orbital parameters and found that tidal heating 
can be strong enough on some planets to trigger a runaway 
greenhouse. In Fig. 2, we show the configurations that pre- 
dict this state around four different hypothetical M dwarfs. 
The IHZ boundaries are the moist greenhouse limits of Selsis 
et al. (2007) and with the same format as Fig. 1. The colored 
curves mark where Fige=Foit Red curves assume the CTL 
model, blue the CPL with discrete rotation states; see Ap- 
pendix E.1. Solid curves assume the Pierrehumbert (2011) 
runaway greenhouse model, and dotted the dry world 
model of Abe et al. (2011); see Fig. Al in Appendix A. For 
these latter worlds, we choose Q=100, as they do not 
have oceans. Thick lines are for a 10 Mgartn planet and thin 
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for 1 Meartn [Abe et al. (2011) only considered a 1 Meartn 
planet]. Note that our choice for the relationship between Q 
and t influences which model predicts more heating. If we 
had chosen the same relationship as that of Heller et al. 
(2011), we would have found that the CTL model predicts 
more heating than CPL; compare their Fig. 5. 

For the lowest-mass M dwarfs, the HZ is significantly re- 
duced due to tidal heating. For masses larger than 0.25 Msun, a 
tidal greenhouse in the IHZ is only possible on large-mass 
planets with very large eccentricities. For lower mass stars, 
our model predicts a tidal greenhouse at low eccentricities 
where tidal theory is most likely valid. Figure 2 shows that a 
planet may have a climate catastrophe due to tide-driven 
overheating even if it is far enough from the star that stellar 
radiative heating alone would not preclude habitability. 

The difference between curves representing equal-mass 
planets is due to the frequency dependence of the CTL 
model. For Earth, the frequency is the mean motion of the 
lunar orbit, but a planet at a=0.035 AU, that is, the middle of 
the HZ, orbits in about 1 week. One could formally adjust 1 
so that it is equivalent to a Q of 10 (see, e.g., Matsumura et al., 
2010; Heller et al., 2011; Eq. E21); then the curves would lie in 
a similar location. The unknown tidal response of terrestrial 
bodies and the absence of an unambiguous translation from 
the CPL to the CTL model motivates our choice of adopting 
present-Earth values for all planets at all frequencies. 

Regions to the left of the colored curves in Fig. 2 can 
produce Tidal Venuses if the planets remain there longer 
than tages. In Fig. 3, we show the evolution of two example 
systems consisting of a 10 Megarth planet orbiting a 0.1 Msun 
star, using both the CPL and CTL models. For the CPL case, 
the planet begins with semimajor axis a=0.04 AU, and e=0.3. 
The planet is assumed to be spin-locked and with zero 
obliquity. This orbit is toward the outer edge of the IHZ and 
with the typical eccentricity of known exoplanets. For the 
CTL case, everything is the same, except the initial a is 0.035 
AU, as the Tidal Venus region lies closer to the star. The top 
panel in Fig. 3 shows the evolution of a, the next panels 
down show e, then insolation Fi,>,1, then tidal heat flux Frage, 
and finally the sum of the insolation and tidal heat flux, Frotai. 
The tidal heat flux generated in the planet in the CPL case 
has decreased to F ..;,=345 W m * at 200 Myr, and the sum of 
tidal heating and insolation reaches F.,i, at 275 Myr. In the 
CTL case, Fridge = Fort at 130 Myr, and Fiotay=Ferit at 275 Myr. 
We conclude that planets orbiting low-mass stars on eccen- 
tric orbits may experience tidal greenhouse conditions long 
enough to become uninhabitable. 

The CPL model predicts a will increase through angular 
momentum transfer with the star, for which we assume a 
rotation period of 30 days, a tidal Q of 10°, a radius deter- 
mined by the Reid and Hawley (2000) relation, and both 
the radius of gyration and Love number of degree 2 are 
0.5. However, this orbital expansion is a result of the discrete 
nature of the CPL model, which only includes four “tidal 
waves,” see Appendix E.1. In Fig. 3, the planet initially spins 
three times per two orbits; hence ¢,,=0, see Eq. E5, elimi- 
nating one of the terms in da/dt. Note that once the eccen- 
tricity drops below \/1/19 the spin rate becomes synchronous 
(see Appendix E.3); then 9 ;=0 but ¢,,=1, and then da/dt <0. 
For this example, that transition occurs at 375 Myr and is not 
shown in Fig. 3. This point illustrates the complexity inherent 
in this commonly used tidal model. 
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FIG. 2. Parameter space of the Tidal Venus. The top left panel is a 0.1 Moy, star, top right 0.15 Msun, bottom left 0.2 Msun, 
and bottom right 0.25 Msun. The grayscale represents the Selsis et al. (2007) IHZ boundaries: from lightest to darkest gray, the 
cloud coverage is 100%, 50%, and 0%, respectively. The colored curves mark where Frige=Ferit. Red curves assume the CTL 
model, blue the CPL. Solid curves assume the Pierrehumbert (2011) runaway greenhouse model and dotted the dry world 
model of Abe et al. (2011). Thick lines are for a 10 Mean planet and thin for 1 Megan. Tidal Venuses lie to the left of these 


curves. 


5. Constraining Observed Planets 


The previous section demonstrated that planets may be- 
come desiccated by tidal heating, but in many cases we will 
be more interested in the possibility that such a condition 
developed on a planet that we have discovered. As tides 
tend to circularize orbits, we may find a planet in the IHZ 


with low eccentricity that experienced a tidal greenhouse 
early on and is hence currently uninhabitable. In this section, 
we describe how to evaluate a known planet’s probability for 
habitability based on past tidal heating. 

To model the tidal heating history of an exoplanet, we 
require knowledge of M,, Rp, 4, e, Mx, R., and age. Exo- 
planets are predominantly discovered by radial velocity and 
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FIG. 3. Evolution of a 10 Megan planet orbiting a 0.1 Ms, star with an initial orbit of a=0.04 AU, e=0.3. Top: semimajor axis 
evolution. Top middle: eccentricity evolution. Middle: insolation evolution. Bottom middle: tidal heat flux evolution. Bottom: total 


surface heat flux (insolation + tidal). 


transit studies that can provide a and e, while the others can 
be modeled from the stellar spectrum. If some of these values 
are unknown, we can often use scaling relations to estimate 
them; see Appendices C—D. From this information, one can 
estimate the tidal heating history of a planet, using the tidal 
evolution equations (Appendix E), and hence estimate the 
probability that the planet has lost its water. As an example, 
consider the hypothetical situation in which a terrestrial- 
scale planet has been discovered in orbit about an M dwarf. 
Further, the planet lies near the inner edge of the IHZ and 
has a low eccentricity. To evaluate past tidal heating, we 
created four systems with different stellar and planetary 
properties, applied reasonable uncertainties to each observ- 


able property, chose non-observable properties via appro- 
priate scaling laws, and modeled the system’s history. In 
particular, we assume that the planet’s mass lies between 1 
and 5 Megarth, has an eccentricity less than 0.1, a tidal Q, (t,) 
in the range 10-100 (64-640), and is tidally locked. We as- 
sume the star’s age lies between 2 and 8 Gyr, and Q, (t.) in 
the range 10° to 10°° (1 to 10°* s) and #g,.=10° years. We 
used the same scaling relations as in Section 4. For the 0.1 
Msun case, a was chosen in the range [0.029, 0.031] AU, for 
0.15 Msyn it was [0.055, 0.065] AU, for 0.2 Ms, it was [0.085, 
0.095] AU, and for 0.25 Mg, it was [0.095, 0.105] AU. For 
each parameter the uncertainty distribution was uniform in 
the quoted range, except for e, which was chosen uniformly 
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in the range —5<logio(e)< —1. We then randomly deter- 
mined the system parameters in these ranges and integrated 
the tidal history backward in time, using the models pre- 
sented in Appendix E, for the randomly chosen age of the 
system. We ignored the possibility of spin-orbit resonance 
capture, which would dramatically alter the history. For each 
stellar mass and tidal model, we simulated 30,000 possible 
configurations. 
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In Fig. 4, we show our results graphically. As we are 
considering a range of masses, we chose to represent the 
runaway greenhouse flux with a 2.5 Mgartn planet; that is, the 
solid curves show where Fiige=Feri¢= 309 W m7. The three 
contours show the probability density for the planet’s loca- 
tion taes= 10° years after the system’s formation. The con- 
tours show where the probability has dropped by 50%, 90%, 
and 99% of the peak value; solid contours are for the CPL 


0.10 
Semi—Major Axis (AU) 


FIG. 4. Orbits of hypothetical planets around M dwarfs after ty... The top left panel is a 0.1 Msy, star, top right 0.15 Msun, 
bottom left 0.2 Msun, and bottom right 0.25 Msun. The grayscale represents the Selsis et al. (2007) IHZ boundaries with the 
same format as Fig. 2. For reference, the red and blue curves show the tidal greenhouse limit for a 2.5 Mgarntn planet with 
Q,=10 (CPL model, blue curve) or t,=640s (CTL model, red curve). Contours denote levels of constant probability density 
(for densities 50%, 90%, and 99% of the peak value) for the initial orbit of the planet: solid corresponds to the CPL model 
(compare to blue curve), dashed to CTL (compare to red curve). In the bottom two panels there has been negligible orbital 


evolution for either tidal model. 
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model, dashed for the CTL. If the probability contours in- 
tersect or even come close to the colored Fyige=F eit Curves, 
then the planet may be a “Habitable Zone Venus,” a planet 
that appears habitable by the IHZ metric but is probably 
more Venus-like than Earth-like. 

The shapes of the regions in Fig. 4 are due to our pa- 
rameter choices and the assumptions implicit in the two tidal 
models. The CPL model does not include eccentricity terms 
to as high an order as the CTL model; hence the CTL model 
predicts more evolution at large e. In the top left panel, this 
effect is seen clearly, as the most likely orbits at 100 Myr in 
the CPL model are at lower semimajor axis than in the CTL 
model. We stress that this example is purely hypothetical, 
and the actual shapes of the contours could be very different 
for actual systems. 

For the 0.1 Msun case (top left), the planet may have spent 
enough time in the tidal greenhouse to be uninhabitable. For 
larger-mass stars, however, the danger of tidal desiccation is 
smaller. In the 0.15 Msun case, the peak in the CPL proba- 
bility density at a=0.08, e=0.45 represents 1% of our simu- 
lations. Although the contours do not cross the CPL runaway 
greenhouse (blue) curve, they do come close; hence there is a 
small chance that our putative candidate is a super-Venus, 
especially if we allow for absorption of stellar radiation. 

However, for M,.>0.15 Mgun, planets with low eccentricity 
probably have always had low eccentricity; that is, the evo- 
lution was negligible. This sharp contrast between 0.1 and 
0.2 Msun occurs because of the steep dependence of tidal 
heating on a. At larger stellar masses, the circular [HZ has 
been pushed away from the reach of fatal tidal heating. 

This experiment is purely illustrative, with arbitrary pa- 
rameters and uncertainty distributions. If a terrestrial planet 
was discovered in the IHZ of a very low-mass star, then this 
methodology could be performed in order to characterize its 
potential to support life. One could derive a value of tues 
tailored to the primary and determine the probability that 
the planet spent more than that in a tidal greenhouse state 
and hence is uninhabitable. In principle, one should also 
allow a range of albedos and include radiation to determine 
the amount of time that the total surface energy flux, Fror= 
Fiiae + Finsor is larger than F,,i:; but we leave such an analysis 
for future work. 


6. Application to GI 667C 


Two or more planets orbit the 0.3 Msun star Gl 667C 
(Bonfils et al., 2011; Anglada-Escudé et al., 2012; Delfosse 
et al., 2012). Planets c and d appear to lie in the IHZ, while a 
third planet, b, lies interior. Planet c is at least 4.5 Mgartn and 
orbits near the inner edge of the IHZ. Planet d, which is 
weakly detected in both studies, lies near the 50% cloud 
cover outer boundary. At 0.23 AU, planet d is too far from 
the star to be subjected to strong tidal heating, at least if it is 
tidally locked. Anglada-Escudeé et al. (2012) and Delfosse et al. 
(2012) proposed different solutions to the system with the 
former setting c’s eccentricity to 0 but stating that it is only 
constrained to be <0.27, while the latter assign its eccen- 
tricity to be 0.34+40.1. The minimum mass estimates are al- 
most identical at 4.25 Meartn (Anglada-Escudé et al., 2012) 
and 4.5 Mrartn (Delfosse et al., 2012). Here, we use the An- 
glada-Escudé et al. (2012) solution, but note that using data 
from (Delfosse et al., 2012) does not change our results. 
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FIG. 5. Comparison of the [HZ to the tidal greenhouse limits 
for Gl 667C c, in a similar format as Fig. 2. Here the thin lines 
correspond to a 4.5 Megan planet and thick to 9 Mgarth. The 
vertical black lines correspond to the uncertainty in eccen- 
tricity of c and d (the latter’s existence remains uncertain). The 
dashed, black curve represents where c’s and b’s orbits cross; 
hence the region to the left is dynamically unstable. 


In Fig. 5, we show the system in the same format as in Fig. 2. 
As this planet was detected via radial velocity data, its true 
mass is unknown. We consider two possibilities here: an 
edge-on geometry in which its actual mass is 4.5 Mean and 
an inclined case in which its actual mass is doubled, 9 Mgartn. 
In Fig. 5, the thin lines correspond to the minimum mass, thick 
to twice-minimum. The vertical extent of the line corresponds 
to the uncertainty in eccentricity. The orbit of c is marked by 
the vertical black line at a=0.123 AU, and d at 0.23 AU. 

The CTL model (red curves) barely intersects the 100% 
cloud cover IHZ at large e. The CPL model penetrates the [HZ 
more significantly but still does not reach the orbit of c. Thus, 
planet c is not currently experiencing a tidal greenhouse, as- 
suming no obliquity and pseudo-synchronous rotation. We 
also see that d is completely safe from the tidal greenhouse. 

But could planet c have been in the tidal greenhouse and 
be desiccated today? The answer is almost assuredly no. In 
this case, we can appeal to the orbital architecture, rather 
than running a suite of Monte Carlo simulations as in Section 
5. The dashed, black line of Fig. 5 shows where the orbit of c 
crosses the orbit of b. Such a configuration is as close as the 
two planets could possibly be and remain stable, assuming 
no mean-motion resonances, which are not detected in this 
system (Anglada-Escudé et al., 2012). This black curve is 
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always exterior to the tidal greenhouse curves out to the 
current orbit of c. Therefore, we may safely conclude that c 
was never in a tidal greenhouse. Furthermore, the planet 
should have become tidally locked within 10’ years (Heller 
et al., 2011), which is much less than tges. Hence, any initial 
burst of tidal heating due to non-equilibrium rotation or 
obliquity was too short to sterilize this planet. Gl 667C c 
remains a habitable planet candidate. 


7. Discussion 


The previous results show where tidal heating can limit 
habitability for planets orbiting low-mass stars. However, 
many phenomena complicate the process and could alter our 
findings. Here, we discuss the results of the previous sections 
in light of our simplifying assumptions, observational re- 
quirements, and theoretical limitations. 

We find that Gl 667C c probably did not lose its water to 
tidal heating, because interactions with other planets prevent 
its eccentricity from being large enough to trigger the tidal 
greenhouse. As more potentially habitable planets are dis- 
covered around low-mass stars, a similar analysis as in 
Section 5 should be undertaken in order to assess the pos- 
sibility that the planet could in fact be dehydrated. As we 
may only be able to spectroscopically characterize a few 
planets with the James Webb Space Telescope (Kaltenegger 
and Traub, 2009; Seager et al., 2009), prioritization of targets 
is crucial, and past and present tidal heating will help de- 
termine the best planet to observe. 

The timescale for e to decay may be smaller than the 
timescale for in situ terrestrial planet formation (Lecar and 
Aarseth, 1986; Wetherill and Stewart, 1989; Lissauer, 1993; 
Kokubo and Ida, 1998; Raymond et al., 2007). On the other 
hand, terrestrial planets could be pushed into such a position 
by a migrating gas giant in about that timescale (Raymond 
et al., 2006; Mandell et al., 2007). It is therefore natural to 
wonder whether such planets can even form, as the tidal 
effects could suppress planet formation or damp out the 
eccentricity of a protoplanet before it is massive enough to 
support clement conditions. Several scenarios suggest that 
fully formed planets can be Tidal Venuses. Orbital in- 
stabilities can excite e. Planet-planet scattering and divergent 
resonance crossing appear to play a role in sculpting many 
planetary systems, including our own (e.g., Weidenschilling 
and Marzari, 1996; Tsiganis et al., 2005; Nesvorny, 2011). 

These phenomena can develop long after planet formation 
has occurred. Recently, it has been suggested that exoplan- 
etary systems in resonance appear systematically younger 
than the general population, suggesting that instabilities can 
even occur after many billions of years (Koriski and Zucker, 
2011). Therefore, we conclude that tidal greenhouses are 
plausible. 

The role of oceans in the tidal dissipation process is clearly 
very important yet also very poorly understood. We have 
used up-to-date information regarding the dissipation of 
energy in Earth’s ocean and assumed that terrestrial exo- 
planets will behave similarly, yet many issues remain out- 
standing. We have ignored the role of inertial waves, which 
provide additional heating of aqueous mantles of icy satel- 
lites (Tyler, 2008, 2011). Therefore, the tidal heating values 
we obtained may, in fact, be too low, further increasing the 
threat to habitability. Regardless, our ranges for Q or 1 
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spanned an order of magnitude, but a wider range of values 
is still possible, as they are complex functions of ocean depth, 
ocean floor topography, and the shape of any continents or 
islands that may be present. As “exobathymetry” seems a 
distant dream, tidal dissipation in Earth-like worlds will re- 
main mysterious for the foreseeable future. 

Further complicating the situation, oceans may evaporate 
long before tages, potentially leading to a cyclical process of 
evaporation and precipitation: After the oceans disappear, Q 
increases and the tidal heat decreases, causing the planet to 
drop out of the tidal greenhouse, so the water rains out, 
reforming the oceans and lowering Q again. Whether such a 
cycle exists is pure speculation, but we note that an analo- 
gous situation can occur with the classic IHZ, where evap- 
orated water forms clouds that increase the albedo, which in 
turn lowers the upward long-wavelength flux from the sur- 
face, and the planet then drops out of the runaway green- 
house. Our choices for Q, t, Feit, and tge, can be revised as 
new observations provide firmer constraints. Of these, taes 
needs the most work, as it is a function of poorly constrained 
host properties and challenging planetary escape processes 
(see, e.g., Tian, 2009). 

A Tidal Venus is an extreme case of tidal heating, over 
2 orders of magnitude more powerful than on Io. We have 
therefore made a considerable extrapolation. Perhaps oceans 
and/or mantles adjust to the increased heating and fail to 
reach that flux. Our choice for Q and t implies that most 
dissipation occurs in the ocean. For the solid interior, 
Béhounkova ef al. (2011) found that tidal heating at about 
Fait leads to a “thermal runaway” for planets that transport 
internal energy through convection. Although they did not 
consider the possibility of advection via volcanism, the 
geophysics of Tidal Venuses requires closer scrutiny. We are 
unaware of any research exploring the physical oceanogra- 
phy of planets undergoing strong tidal heating, suggesting it 
is an interesting topic for future research. 

Our treatment ignored mutual gravitational interactions 
between planets or large satellites, which can significantly 
alter the evolution. Gravitational perturbations can pump up 
eccentricities and obliquities to nonzero values and may be 
able to modify the spin period. From Fig. 2, we see that 
planets near the inner edge of the HZ may be in a tidal 
greenhouse for eccentricities less than 0.02. Other bodies in 
the system can easily perturb eccentricities to larger values; 
hence planets in multiple systems may be especially sus- 
ceptible to a tidal greenhouse, as this driven eccentricity can 
be maintained for arbitrarily long timescales. As individual 
systems are discovered, this point should be revisited. 

On the other hand, orbital stability arguments could allow 
us to preclude the tidal greenhouse in multiplanet systems, 
as for Gl 667C. When the eccentricity required for the tidal 
greenhouse is so large that the system is unstable, then one 
can safely exclude the Tidal Venus state. Hence, the presence 
of additional companions provides critical information when 
assessing habitability. 

Many other critical phenomena were also left out, such as 
atmospheric erosion by flaring, stellar activity, and magnetic 
dynamo generation (Khodachenko et al., 2007; Lammer et al., 
2007, 2010; Tian, 2009; Segura et al., 2010; Driscoll and Olson, 
2011). Tidal heating provides an interesting counterbalance 
to atmospheric erosion, as it may increase the outgassing 
rates and maintain a permanent atmosphere. The outgassing 
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and escape need to remain in a balance, as wild pressure and 
density fluctuations will undoubtedly alter the biosphere; 
but in principle, tidally driven outgassing could reduce the 
danger of atmospheric removal. On the other hand, more 
intense outgassing without commensurate drawdown by 
processes like the carbonate-silicate cycle could increase the 
threat of a moist greenhouse. Magnetic fields may slow 
down atmospheric loss, but tidal heating may decrease the 
dynamo. If magnetic fields are generated by convection in the 
outer core between a hot inner core and a cool mantle (Olson 
and Christensen, 2006; Driscoll and Olson, 2011), then tidal 
heating of the mantle may suppress magnetic field generation 
(Stevenson, 2010). This issue has not yet been explored for 
tidally heated exoplanets, but it could play a major role in 
habitability. Future work should couple outgassing and es- 
cape rates in order to determine how the two interact. 


8. Conclusions 


We have shown that tidal heating of some exoplanets may 
exceed the threshold of the runaway greenhouse, the tradi- 
tional inner edge of the IHZ. We find that for stars with 
masses < 0.3 Msun, planets in their IHZs with low eccentricity 
can be uninhabitable regardless of insolation. We have thus 
fundamentally revised the HZ boundaries for planets on 
eccentric orbits. Unlike insolation from main sequence stars, 
tidal heating at a desiccating greenhouse level may drop off 
rapidly, but not so rapidly as to preclude the possibility that 
a planet’s entire inventory of water can be lost permanently; 
see Fig. 3. These planets will be uninhabitable regardless of 
future tidal heating, that is, a planet found with minimal 
tidal heating today may still have experienced sufficient 
heating for sufficient duration to render it uninhabitable. 
Additional planetary companions are important: They can 
drive eccentricity and sustain a tidal greenhouse, or they can 
be used with stability arguments to rule out an early tidal 
greenhouse. 

Traditionally, habitability models have focused on insola- 
tion, implying that the star is the most important aspect of 
habitability. We have shown that in some circumstances tidal 
effects are more important in determining the inner edge of the 
HZ. Planetary habitability is a function of the star, the planet, 
and the planetary system. The Kasting et al. (1993) THZ has 
served well as a guide but is insufficient. Combining all the 
processes relevant to habitability into a single model is a 
daunting challenge to say the least, but a proper assessment of 
a planet’s potential for habitability relies on a wide diversity of 
properties, some of which will not be observable any time 
soon. Nevertheless, the prospect of identifying an inhabited 
planet is strong motivation. Moreover, the high cost in time, 
money, and resources required to establish a planet as in- 
habited demand that we use these resources efficiently. In this 
study, we have compiled numerous tidal processes and em- 
pirical relationships so that at least the tidal effects predicted 
by linear theory may be applied to terrestrial planets in the 
IHZs of low-luminosity hosts. 


Appendix A. Relations for the Inner Edge 
of the Habitable Zone 


The critical flux for the runaway greenhouse depends 
weakly on the acceleration due to gravity, g, in the planetary 
atmosphere. Flux is emitted to space from the level of optical 
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depth unity. The optical depth of a layer of absorbing gas 
depends to first order on the mass of the layer, so with higher 
g the pressure at the base of a constant mass layer will be 
higher. For the Simpson-Nakajima limit, the T—p structure of 
the atmosphere is the saturation vapor pressure curve, so 
higher p implies higher T and permits a higher critical flux. 
This positive relationship is offset to some extent by pres- 
sure-broadening of absorption, so at higher g (and higher p), 
less mass of water is required to provide an optical depth of 
unity (Pierrehumbert, 2010; Goldblatt and Watson, 2012). 
Pierrehumbert (2010) derived a semi-analytical relation- 
ship for the critical flux due to the Simpson-Nakajima limit: 


4 


l 


Rin(P. a] 


where | is the latent heat capacity of water, R is the universal 
gas constant, Po is the pressure at which the line strengths are 
evaluated, g is the gravitational acceleration in the planetary 
atmosphere. x is a gray absorption coefficient, and A is a 
constant of order unity; Pierrehumbert (2010) fits these to 
numerical runs of a radiative transfer code with «=0.055 and 
A=0.7344. P, is a scaled pressure given by 





Feat =Ao (Al) 


l 
P= Pee (A2) 
where Tyep=273.13 K and P,.¢=610.616 Pa are points on the 
saturation vapor pressure curve of water, taken as the triple 
point. This calculated limit is shown in Fig. A1. 
For convenience, we use an empirical relationship for the 
moist greenhouse limit derived by Selsis et al. (2007): 
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The outer edge (not due to a desiccating greenhouse) can be 
expressed in a similar format as 
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FIG. Al. Surface flux required to trigger a runaway 


greenhouse (F.,i) aS a function of planet mass. The solid line 
is the relation for wet planets from Pierrehumbert (2011), and 
the cross is the limit from Abe et al. (2011) for a dry planet. 
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(A4) 
In these equations, J, and /,u4 are the inner and outer edges 
of the IHZ, respectively, in astronomical units; linsun and 
loutSun are the inner and outer edges of the IHZ in the solar 
system, respectively, in astronomical units; ain =2.7619 x 107° 
AU/K, bi, =3.8095 X 107? AU/K?, a5,4= 1.3786 107 * AU/K, 
and bou.=1.4286x10~”? AU/K? are empirically determined 
constants; and L, and Ls,, are the stellar and solar lumi- 
nosity, respectively. T, =Tegg — 5700 K, where Ter is the “ef- 
fective temperature” of the star 


a L, 1/4 
ot \ 4nR? 


where R,, is the stellar radius and oa is the Stefan-Boltzmann 
constant. Usually, however, the IHZ is couched in terms of 
the stellar mass M.,, and a; hence relationships between L,, 
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FIG. A2. IHZ boundaries for different combinations of 
mass-radius and mass-luminosity relations and using the 
three different cloud cover IHZ limits from Selsis et al. (2007). 
Except for the inner edge near M,.=0.2 Moun, the effective 
temperature does not affect the boundaries (note the nu- 
merous curves that are visible). Six combinations are plotted, 
but most are invisible as they are on top of each other. At low 
M,.,, the more interior [HZ limits assume the RHOO mass- 
luminosity relationship; see Appendix C. 
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M,, and R, are needed to produce such an expression. Em- 
pirical and theoretical relations between these properties are 
reviewed in Appendix C. 

The values of linsun and [outsun are therefore the key 
parameters in the identification of the edges of the IHZ. We 
consider three criteria identified by Selsis et al. (2007): (1) 
0% cloud cover, (2) 50% cloud cover, and (3) 100% cloud 
cover. Selsis et al. (2007) give values of linsun of 0.89, 0.72, 
and ~0.49 AU for the three possibilities, respectively. For 
the outer edge, Selsis et al. (2007) gave values of Ijutsun Of 
1.67, 1.95, and 2.4 AU for the three cloud cover models, 
respectively. We arbitrarily choose the 50% cloud cover 
case to be the limits of the IHZ. This choice is the middle of 
the possibilities and roughly corresponds to the cloud 
cover on Earth. 

However, as Selsis et al. (2007) noted, 


the effect of the spectral type on the albedo, included in 
Eqs. (3-4) as a quadratic function of (Teg — 5700), was 
estimated only for a cloud-free atmosphere. Since the 
reflectivity of clouds is less sensitive to wavelength, 
this quadratic term may not be valid to scale the 
boundaries of the HZ for planets covered by clouds. 


Furthermore, the assumption that Selsis et al. (2007) made, 
that clouds do not affect the IR flux from a planet at the inner 
edge of the HZ, is not strong, so the cloudy results should be 
interpreted with caution. Furthermore, assumptions used in 
the models from which Selsis et al. (2007) derived their re- 
sults assume fast planetary rotation, so it will not be a good 
approximation for most tidally locked exoplanets. 

Figure A2 shows the range of the IHZ for the M dwarf 
mass range for the empirical relations presented in Appendix 
C. For very low-mass stars (<0.1 Msun), the position of the 
THZ may differ by 50% depending on the relationships in- 
voked. However, for larger stars, the [HZ is mostly inde- 
pendent of the chosen mass-luminosity and mass-radius 
relationships. 


Appendix B. The Desiccation Timescale 


As discussed in Section 2, the details of mass loss on a 
planet near the inner edge of the HZ of a very low-mass star 
are largely unknown. Therefore, the timescale for desiccation 
is difficult to constrain. In this appendix, we describe a 
simple model for tye, and justify the value used in the main 
text of 10° years. 

We use the atmospheric mass loss model described by 
Erkaev et al. (2007), which improves the standard model by 
Watson et al. (1981). In the Watson et al. picture, a layer in 
the atmosphere exists where absorption of high-energy 
photons heats the particles to a temperature that permits 
Jeans escape. Photons in the X-ray and XUV have the en- 
ergy to drive this escape on most self-consistent atmo- 
spheres. In essence, the particles carry away the excess solar 
energy. Watson et al. (1981) calculated that Venus would 
lose its water in 280 Myr. The actual value is probably less 
than that, as they were unaware that XUV emission is larger 
for younger stars, as the spacecraft capable of such obser- 
vations had yet to be launched. Furthermore, their model 
did not consider the possibility of mass loss through La- 
erange points, or “Roche lobe overflow.” Erkaev et al. 
(2007) provided a slight modification to the Watson et al. 
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(1981) model that accounts for this phenomenon for hot 
Jupiters. The physics of mass loss from a terrestrial atmo- 
sphere may be different than that on a giant, but as we are 
interested in objects that are experiencing tidal forces, we 
use the Erkaev et al. (2007) model, 


dM, 
dt 


R2R,eF 
oa ee (B1) 


where R, is the radius of the atmosphere at which the op- 
tical depth for stellar XUV photons is unity, ¢ is the effi- 
ciency of converting these photons into the kinetic energy of 
escaping particles, Fyyy is the incident flux of the photons. 
The parameter kiige represents the reduction in energy re- 
quired for escape of a particle due to the star’s gravitational 
pull: 
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is the ratio of the “Hill radius” to the radius at the absorbing 
layer. From Eq. B1 it is trivial to show that 
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where mr is the total mass of all the hydrogen atoms that 
must be lost. 

The values of Fxyy and ¢ are therefore the key parame- 
ters, as they represent the magnitude and efficiency of the 
process. Unfortunately, very little is known about the for- 
mer for very low-mass M dwarfs with ages of a few hun- 
dred million years. The stars are very faint and, moreover, 
prone to flaring. However, there is a secular decrease 
in activity as the stars age (West et al., 2008), but the data 
are still too sparse to allow a satisfactory fit. Recent 
work connected X-ray sources from space telescopes (e.g., 
GALEX) and ground-based surveys (e.g., Sloan Digital Sky 
Survey) and found that older low-mass stars tend to emit 
about 10° to 10° of their energy in the XUV in quiescence 
(West, personal communication). During flares, that emis- 
sion can increase by an order of magnitude (Hawley and 
Pettersen, 1991; Kowalski et al., 2012). As all stars are more 
active when they are young, we should expect the XUV 
emission to be substantially larger at those times. 

The value of ¢ must be less than 1, as not all the absorbed 
energy removes particles. Observations of hot Jupiters sug- 
gest values of ¢ of 0.4 (Yelle, 2004; Lammer et al., 2009; 
Jackson et al., 2010), while on Venus e~0.15 (Chassefiére, 
1996). However, our model also requires these photons to 
dissociate the water molecules, hence we should expect ¢ on 
a Tidal Venus to be much less than on a hot Jupiter with a 
predominantly hydrogen atmosphere. For reference, using 
the assumptions of Watson et al. (1981), we find ¢=1.7x 10° * 
implies tges=280 Myr for Venus. 

In Fig. B1, we show the desiccation timescale for an Earth- 
like planet orbiting a 0.1 Msun star at 0.03 AU for a range of 
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biting a 0.1 Msun star as a function of the fraction of the lu- 
minosity in the XUV and the efficiency of converting that 
energy into escaping particles. Contour lines are logio(taes/Yr). 
The planet initially has the same water mass fraction as the 
Earth and a semimajor axis of 0.03 AU. 


XUV luminosities and «¢. The initial mass of the planet in 
hydrogen is 1.4 x 107° Meartn, the current value of Earth and 
the same value used by Watson et al. (1981), and R, =1.02 R,,. 
The luminosity of the star is 0.0008 Lsu,, and the planet lies 
inside the 0% cloud cover HZ (see Appendix A). The contour 
lines show logio(tges). For e~10-*, a few percent of the lu- 
minosity would need to be in the XUV for tages to be 10° years, 
while hot Jupiter escape rates would only require>10~° of 
the luminosity to be in the XUV. 

The environment near a forming star is very energetic; 
therefore we should expect other phenomena, such as abla- 
tion and magnetic fields, also erode the atmosphere. There- 
fore, the actual values of tye, are assuredly shorter than those 
presented in Fig. Bl. We conclude that tg.,=100 Myr is a 
reasonable choice, while explicitly acknowledging a large 
range is possible. Finally, we note that simulations should 
couple the evolution of the XUV output to the tidal orbital 
evolution of the planet in order to calculate tg., for a given 
system. Jackson et al. (2010) created such a model for 
CoRoT-7 b (Léger et al., 2009), which orbits close to a solar- 
mass star. For such stars the generic evolution of their XUV 
evolution is known (Ribas ef al., 2005); hence such an 
analysis was quite revealing. Given the poor constraints on 
XUV evolution for M dwarfs, we eschewed such coupling, 
deferring that analysis to future work. 
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Appendix C. Masses, Radii, and Luminosities 
of M Dwarfs 


In this appendix, we review relationships between mass, 
radius, luminosity, and effective temperature for M dwarfs 
(M,.<0.6 Msyn), and their effects on the limits of the IHZ. 
Measurements of these fundamental stellar properties are 
challenging, and several studies have produced empirical 
and theoretical relations. 

We first consider the mass-radius relation. We include 
empirical models by Gorda and Svechnikov (1999) (GS99), 
Reid and Hawley (2000) (RH00), Bayless and Orosz (2006) 
(BO06), as well as the theoretical models of Baraffe et al. 
(1998) (B98). This latter model is parameterized in terms of 
the metallicity as defined by astronomers: 


[Fe/H] = log,, [Fe/H], — log,, [Fe/H]sun (C1) 


in other words, it is the difference between the ratio of iron to 
hydrogen in a star to that of the Sun. B98 consider solar 
metallicity stars ([Fe/H]=0) and subsolar ([Fe/H]=-—0.5). 
GS99 and BO06 present analytic fits that are shown in Table 
C1. For RHOO, we take their data in Table 4.1 and derive a 
third-order polynomial fit, y=bo+b,x+box*+b3x°, using a 
Levenberg-Marquardt minimization scheme, and present the 
results in Table Cl. The top two panels of Fig. C1 show the 
relationships graphically. 

Next, we examine the mass-luminosity relationships de- 
tailed in RHO0, Scalo et al. (2007) (SO7) and B98. For the 
former, we again fit their data with a third-order polynomial. 
The fits are listed in Table Cl and shown graphically in the 
bottom panels of Fig. C1. 


Appendix D. Radii of Terrestrial Planets 


In this appendix, we review published studies that provide 
analytic formulae relating terrestrial planet mass and radius. 
Throughout the study, we appealed to the relationship de- 
rived by Sotin et al. (2007), who considered two types of 
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planets: Earth-like composition and “ocean planets” that are 
50% HO by weight. The scaling relationship for the former is 
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and for the latter 
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Fortney et al. (2007) parameterized the radii (M, >0.01 Meartn) 
in terms of the ratio of ice to rock, fice, or rock to iron, frock. They 
assumed that these two combinations are the most likely for 
terrestrial exoplanets and derived the following relationships: 
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+ (0.2337 frock + 0.4938) log, 





Earth 
+ 0.3102 frock + 0.7932 


In these equations, fi-e=1 corresponds to a pure water, fice=0 
and frock=1 to a pure rock (silicate) planet, and fioxK=0 is a 
planet made of pure iron. Note that Eqs. D3—-D4 are analytic 
fits to complex models and, hence, predict slightly different 
values for R, for a planet made of rock (fice=0 and frock = 1). 


TABLE Cl. SCALING RELATIONS FOR M DWARES 








Xx V bo b, b> b3 y° Ref. 
logio(M.. / Msun) logio(R. i RGun) 0.1 1:03 0 0 n/a GS99 
logio(M./Msun) logio(R«/Rsun) 0.0128 2.185 3.1349 1.903 0.0136 RHO0 

M,./Meun R./Rsun 0.0324 0.9343 0.0374 0 n/a BO06 
logio(Mx/Msun) logio(R«/Rsun) 0.0372 1.5 1.09 0.5361 1.72x10~* B98([Fe/H]=0) 
logio(Mx/Msun) logio(R«/Rsun) 0.0664 1.685 1.398 0.7139 4.38x10~* B98([Fe/H]= -— 0.5) 
logio(Re/Rsun)  logio(Mx/Msun) —0.0971 0.971 ~2.501x107°  -1.34x107° 1.2x107' Gs99 
logio(R«/Rsun) logio(Mx/Msun) —— 0.1424 1.568 ~ (0.2342 ~0.5581 0.0331 RHO00 

Roe... M,./Msun ~0.03477. ‘1.071 =8171X10 = =—0.0412 2.83x10°? BO06 
logio(Re/Rsun)  logio(M./Msun) —0.03 0.507 ~1.156 ~ 0.5978 8.88x10~* B98([Fe/H]=0) 
logio(Re/Rsun)  logio(M./Msun) —0.0406 0.4537. -1.211 ~ 0.6427 0.00139 B98([Fe/H] = —0.5) 
logio(M./Msun) logio(Ls/Lsun) 0.2984 8.7116 11.562 6.241 0.0155 RHO0 

logio(M.. / Msun) logio(L. Figen) 0.065 7.108 8.162 4.101 n/a S07 
logio(M./Msun) logio(Ls/Lsun) 0.0182 7.228 9.357 5.281 0.039 B98([Fe /H] =0) 
logio(M./Msun) logio(Ls/Lsun) 0.4377 8.908 12.21 6.896 0.094 B98([Fe/H] = —0.5) 
logio(Le/Lgun) — logio(M/Msun) —0.3076 —-0.5 — 0.4504 — 0.06852 0.00287 RHO00 
logio(Le/Lsun) — logio(M./Msun) —0.3536 -0.5463  -0.4422 ~ 0.06243 0.224 S07 

logio(Le/Lsun) — logio(Mx/Msun) —0.04 =A ~ 0.2765 ~0.045 0.039 B98([Fe/H] =0) 
logio(Le/Lsun)  logio(Mx/Msun) —0.03 0.00509  -0.2135 ~ 0.0368 0.0145 B98([Fe/H] = —0.5) 





n/a means y* values are not available. Only reprinted are fits that were derived in the stated work. 
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FIG. C1. Scaling relations for M dwarfs. Top: Mass-radius relations. Black curves are empirical relations: solid from BO06, 
dashed from RHO0, and dotted from GS99. Red curves are theoretical curves from B98: solid assumes [Fe/H] =0, dashed [Fe/ 
H]= —0.5. Bottom: Mass-luminosity relations. Black curves are empirical relations: solid from RH00 and dotted from S07. Red 


curves are theoretical curves from B98: solid assumes [Fe/H]=0, dashed [Fe/H] = —0.5. 


Earth has frock=0.67; hence Eq. D4 becomes R,/Rearth= Rp (1 +£0.56 fee) My 0.262(1 — 0.138 fice) (5) 
0.1372[logio(Mp/Meartn) | + 0.6504 logio(Mp/Megartn) + 1.0010. Rea ee auas |) oe 

Valencia et al. (2007) considered planets in the range 1-10 
Mearth and parameterized the radius in terms of fice (see also 
Valencia et al., 2010). They found 





Seager et al. (2007) considered a wide range of planetary 
compositions and developed a general, but complicated, 
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mass-radius relationship. Although they developed several 
formulae, here we only present one. They found that by 
casting M, and R, in terms of “scaling values,” all planets 
follow the relationship 

logig ts =n +k logy ms — kom’ (D6) 
where k, = — 0.20945, k,=0.0804, and k3=0.394. The param- 
eters m, and r, are the scaled mass and radii, respectively, 
and are defined as ms=m/m, and r;=r/1r;, where the de- 
nominators are a scaling factor that depends on composition 
and thermal effects. We will consider four compositions. For 
pure iron planets, 1, =5.8 Mgartn and 1, =2.52 Reartn. For pure 
rock (perovskite) planets, m,;=10.55 Mgarth and 1; =3.9 Rearth. 
For pure water planets, m,=5.52 Mearth and 7; =4.43 Regarth. 
Finally, for an Earth-like planet (30% iron, 70% perovskite), 
m,=6.41 Meaty and r,=2.84 Rpatn. Note that this formula- 
tion is only valid for m,<4. 

The last relationship we review is that of Grasset et al. 
(2009), who cast their model in terms only of fice. Their for- 
mulation is complicated, and rather than reproduce it here, 
the reader is referred to their paper. We do not include the 
recent models put forth by Swift et al. (2012), as they did not 
include formulae, but they did find that their results are 
consistent with those of Fortney et al. (2007) and Valencia 
et al. (2010) for rock and iron planets. 

In Fig. D1, we compare these relationships for 0.3-30 
Mearth and for a variety of compositions. Not all the rela- 
tionships span that mass range, nor do they encompass all 
compositions. Furthermore, although we group planets into 
four types, not all types are identical. For example, we labeled 
the Sotin et al. (2007) ocean planet as a “Water” planet, but it 
is only 50% water, while the other three are 100% water. 
Nonetheless, the plot shows that if only mass is known, the 
radius of a planet can only be known to a factor of 22.5 (of 
course, its terrestrial nature is also unlikely to be known). As 
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FIG. D1. Mass-radius relationships for terrestrial planets of 
varying compositions. Solid curves assume a roughly Earth- 
like composition, dashed are all silicates, dotted are at least 
half water, and dot-dashed are pure iron. Black correspond 
to Fortney et al. (2007), red to Sotin et al. (2007), blue to 
Grasset et al. (2009), green to Valencia et al. (2007), and or- 
ange to Seager et al. (2007). 
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tidal heating scales as R°, this translates to a difference in 
heating of 30-100, assuming identical Q or t. Note that the 
solid red curve [Earth-like planet from Sotin et al. (2007)] has a 
typical shape for its class; hence the results presented here do 
not vary much if a different model is chosen, unless a sig- 
nificantly different composition is invoked. 


Appendix E. Tidal Theory 


For our calculations of tidal evolution, we employ “equi- 
librium tide” models, originally conceived by Darwin (1880). 
This model assumes the gravitational potential of the tide 
raiser can be expressed as the sum of Legendre polynomials 
(1.e., surface waves) and that the elongated equilibrium shape 
of the perturbed body is slightly misaligned with respect to 
the line which connects the two centers of mass. This mis- 
alignment is due to dissipative processes within the deformed 
body and leads to a secular evolution of the orbit as well as 
the spin angular momenta of the two bodies. As described 
below, this approach leads to a set of six coupled, nonlinear 
differential equations, but note that the model is, in fact, linear 
in the sense that there is no coupling between the surface 
waves which sum to the equilibrium shape. A substantial 
body of research is devoted to tidal theory (e.g., Hut, 1981; 
Ferraz-Mello et al., 2008; Wisdom, 2008; Efroimsky and Wil- 
liams, 2009; Leconte et al., 2010), and the reader is referred to 
these studies for a more complete description of the deriva- 
tions and nuances of tidal theory. For this investigation, we 
will use the models and nomenclature of Heller et al. (2011), 
which are presented below. 


E.1. The constant phase lag model 


In the CPL model of tidal evolution, the angle between the 
line connecting the centers of mass and the tidal bulge is 
assumed to be constant. This approach has the advantage of 
being analogous to a damped driven harmonic oscillator, a 
well-studied system, and is quite commonly utilized in 
planetary studies (e.g., Goldreich and Soter, 1966; Greenberg, 
2009). In this case, the evolution is described by the following 
equations: 
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where e is eccentricity, t is time, a is semimajor axis, G is 
Newton's gravitational constant, M,; and Mp are the two 
masses, R,; and R> are the two radii, w; are the rotational 
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frequencies, w; are the obliquities, and n is the mean motion. 
The quantity Z’ is 


R? 1 
vas — 3G7ky iMe(Mi + M;) re (E5) 


where kz; are the Love numbers of order 2 and Q; are the 
“tidal quality factors.” The parameter ¢; is 


oo as 
= rR, oan 
‘GM; 





(E6) 


where i and j refer to the two bodies and r, is the “radius of 
gyration”; that is, the moment of inertia is M(r, R)*. The signs 
of the phase lags are 


£0,i = 220: a 2n) 
EL i= 2d (20; —3n 
0 220) = n) 


i= 2(n) (E7) 
63,4 — 2d (a; — 2n) 
é9,; = X (aj) 


with 2(x) the sign of any physical quantity x, thus 2(x)= 
ae Ly tO, 

The tidal heating of the i" body is due to the transfor- 
mation of rotational and/or orbital energy into frictional 
heating. The heating from the orbit is 


, 7, 147 
Fob i= — Xx (4e0,i + e7[ — 20e0,; + — 1; 
8 2 (E8) 


| 
+ 58,1 — 3e5,] — 4 sin’ (W,)[¢0, i — 8,1) 


and that from the rotation is 


27,0; 
Erot,i = — —— X (460, + &7[ —2080, i + 4981, i + £2, 
ot, 3 7 (4é9,; + e*[ —2080,; + 4981, + €2,1] (E9) 
+2 sin’ (W;)[ —2e0,i + 8,1 + £9, i) 
The total heat in the i‘ body is therefore 
*~CPL . . 
Eide,i — sire a Evo) > 0 (E10) 


The rate of evolution and amount of heating are set by three 
free parameters: Q, ky, and r,. We discuss the ranges of these 
parameters below in Appendix E.4. 

Goldreich (1966) suggested that the equilibrium rotation 
period for both bodies is 


CPL _ - 


a Bit 
= 1+9.5e? ee) 


Murray and Dermott (1999) presented a derivation of this 
expression, which assumes the rotation rate may take a 
continuum of values. However, the CPL model described 
above only permits four “tidal waves” and, hence, does not 
permit this continuum. In Fig. El, we compare Eq. E11 to 
the equilibrium values predicted by Eqs. E1—-E4 as a func- 
tion of eccentricity for Gl 581 d, that is, an orbital period of 
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FIG. El. Comparison of the equilibrium spin periods in the 


CPL framework. The solid line is the value predicted by our 
model, Eqs. E1—-E4, the dashed is the “continuous” model of 
Goldreich (1966). 


66 days, a mass of 5.6 Mgartn, and a radius of 1.6 Rearth 
(Mayor et al., 2009). As expected, E7 predicts a continuous 
range of periods, whereas our model (solid curve) does not. 
The jump at e=,/1/19 occurs at the w/n=1.5 transition 
and nearly reaches the value predicted by Eq. E11. The next 
phase jump at e= ,/1/19 ® 0.32 is not included in the CPL 
model; hence it breaks down at larger values. However, at 
larger e values, coupling between waves and other nonlinear 
effects may dampen the role of tides. Therefore, the evolution 
at larger e predicted by the CPL model may not be qualita- 
tively correct. We urge caution when interpreting CPL results 
above e=0.32. 

Thus, the CPL model can be considered from two per- 
spectives. When modeling the evolution of a system, one 
should use the discrete spin values for self-consistency, that 
is, aS an initial condition, or if forcing the spin to remain 
tidally locked. However, if calculating the equilibrium spin 
period separately, the continuous value of Eq. E7 should be 
used. We refer to these rotations as “discrete states” and the 
“continuous state.” Note that this point marks a difference 
between our implementation of the CPL model and that 
described by Heller et al. (2011). 


E.2. The constant time lag model 


The CTL model assumes that the time interval between 
the passage of the perturber and the tidal bulge is constant. 
This assumption allows the tidal response to be continuous 
over a wide range of frequencies, unlike the CPL model. 
But, if the phase lag is a function of the forcing frequency, 
then the system is no longer analogous to a damped driven 
harmonic oscillator. Therefore, this model should only be 
used over a narrow range of frequencies (see Greenberg, 
2009). However, this model’s use is widespread, especially 
at high e, so we use it to evaluate tidal effects as well. This 
model predicts larger tidal heating and evolution rates at 
high e as any coupling between tidal waves is ignored. 
Therefore, the CPL and CTL models probably bracket the 
actual evolution. 


TIDAL VENUSES 


The evolution is described by the following equations: 


























de 11lae fale) a; 18 fa(e) 
di 2GMiMb sc 2 Zoo Ye) 1%) 
da Ww 2 fle ao file) 
dt GMiMp ON (¢ 08 (Yi) ai pre) n ) —) 
dy; = Z; sin (W;) — Gil fale) @ — , fale) 
dt ~ 2Mir2 ;R?no ({cos Ww) 1 Ben 2) 
(E14) 








r= aap (208i) Sa — [1+ cos? yy] OS) 








dt = 2Mir, jRin pe) pr(e) 0 
(E15) 
where 
R° 
Zi = 3G7ky, iMe(Mi + Mj), Ti (E16) 
and 
pe)=v1-8 
OL». 209 A 185 06 25 :: 
A@=1+ yet set eet ge 
ills 45 5 
af cad = 
fo(e) = Po Pe ae 
15 15 5 
fa(e)=1+ Sores ae (E17) 
e) il 
fale) =1+ se + se 
fs(e)=1 43¢ + =o! 
The tidal heating of the i body is therefore 
“CTI, file) fro(e) Oj 
EepeeZ ae ) ~ 22 (¢ ) cos (Yj) — 
(E18) 


1+ cos? (W;)| fs) (a; 
+| 2 Fala) 


Note that Eqs. E6 and E14 include terms due to both rotation 
and obliquity, which were not discussed in the main text. 
Should information regarding these parameters become 
available for particular planets, they should be included in 
estimates of tidal heating and habitability. 

It can also be shown that the equilibrium rotation period 
of each body is 


1+ cos’ wy) 
ie pe fs(e)( El 
C= cr (E19) 
which for low e and w=0 reduces to 
Pp 
CTL. — 
1 ~ Tn6e2 (E20) 
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There is no general conversion between Q,, and t,. Only if 
e=0 and w,=0, when merely a single tidal lag angle ¢,, exists, 
then 

Qp © 1/(2|n — @p|tp) (E21) 
as long as 1 — Wp, remains unchanged. Hence, a dissipation 
value for an Earth-like planet of Q,=100 is not necessarily 
equivalent to a tidal time lag of 638s, so the results for the 
tidal evolution will intrinsically differ between the CPL and 
the CTL model. However, both choices are common for the 
respective model. 


E.3. Numerical methods 


The rates of change of the six parameters affected by tidal 
theory are all different. The orbital parameters tend to evolve 
several orders of magnitude more slowly than the planetary 
spin parameters, at least for the cases we are interested in. 
Therefore, we implement a dynamical time-step scheme 
in our integration of tidal evolution. At the beginning of 
each time step, we calculate the derivatives and set the time 
step to 


(E22) 


where 7 is a constant, which we set to 0.01. We find that 
larger values of 7, even up to unity, produce results similar 
with 7=0.01, but at our more conservative choice, we find 
almost no difference with even smaller value of y; that is, 
the solution has converged. This implementation increases 
the speed of the code by several orders of magnitude in 
some cases. For example, if considering a planet with spin 
properties similar to Earth at 0.01 AU, the timescale to 
damp out the obliquity is of order years, but the orbit 
evolves in millions of years (depending on e). Integrating 
this system for billions of years with a fixed time step 
chosen to resolve the obliquity evolution will take many 
orders of magnitude more time than using this scheme. 
Dynamical time stepping is especially important for planets 
that merge with their host star (Jackson et al., 2009; Levrard 
et al., 2009). 

In many cases, e and W can damp to arbitrarily small 
values; hence de/dt and dy/dt take on meaningless values. In 
some cases, their respective timescales as defined above can 
become very short, even though they are effectively 0. This 
situation destroys the speed advantage gained by utilizing 
E18. We therefore set a floor for these values at 107 '°. When 
they reach this level, we set them to 0 and do not include 
them in the determination of At. 

Finally, care must be taken as a planet approaches tidal 
locking. In many cases, an integration will bounce back and 
forth above and below the equilibrium period. This possi- 
bility is a natural result of the discrete nature of numerical 
integration but, again, can diminish the benefits of a dy- 
namical time step. Therefore, if a planet comes within 10% of 
the equilibrium period, we set the spin period to the equi- 
librium value. If the orbit is eccentric and evolving, after tidal 
locking, we force the planet to remain locked; that is, the spin 
period obeys Eq. E11 or E20. 
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E.4. Tidal response in celestial bodies 


Exoplanets are expected to have an even wider diversity 
of compositions and structures than those in our solar system 
(see Léger et al., 2004; Raymond et al., 2004; O’Brien et al., 
2006; Bond et al., 2010); thus their response to tidal processes 
probably varies greatly. We therefore choose to keep our 
model simple and set ky=r,=0.5 for all planets. 

Earth is the only body in the Solar System with a signifi- 
cant amount of water on its surface and is also the planet 
with the most detailed picture of its tidal processes. The 
majority of tidal dissipation on Earth occurs in the oceans, 
and we assume the same will occur on terrestrial exoplanets. 
Tidal dissipation in Earth’s ocean was directly measured 
from TOPEX and Poseidon satellite data by Egbert and Ray 
(2000), who found that 50-75% of dissipation occurs in and 
around shallow seas, with the remainder in the open ocean, 
presumably due to turbulence caused by topography. Recent 
computer simulations have successfully reproduced this 
observation (Jayne, 2001; Nycander, 2005; Kelly et al., 2010). 
For a more complete treatment of this process, consult the 
work of Garrett and Kunze (2007). Therefore, dissipation in 
planetary oceans is a complex function of continental shapes, 
seafloor topography, and the forcing frequency. 

Lunar Laser Ranging, in which the timing of laser pulses 
bounced off the Moon is measured, has provided very pre- 
cise measurements of the rate of the Moon’s recession from 
Earth (Dickey et al., 1994). If this recession is due only to tidal 
dissipation within Earth, then Earth’s tidal quality factor is 
Orarth=12 (Yoder, 1995). This value is lower than most 
measurements in the Solar System, which are more uncertain 
(see below). However, given the presence of surface oceans, 
this extra dissipation should be expected. On the other hand, 
extrapolating backward from the current orbital expansion 
rate and employing the CPL model, one finds the Moon was 
at Earth’s surface ~2 Gyr ago (MacDonald, 1964). The 
moon-forming impact has been dated to about 4.5 Gyr ago; 
hence Earth’s current Q appears too low. This discrepancy 
may be explained by four possibilities: (1) The current Q is 
anomalously low, and in the past different continental con- 
figurations and/or ocean floor topography did not permit 
such a high dissipation; (2) The tidal effects are a function of 
frequency, and as Earth’s rotation has slowed, the dissipation 
rate has increased; (3) Perturbations from other planets af- 
fected the evolution (Cuk, 2007); or (4) Linear theory is over- 
simplified. Probably all these possibilities play a role and make 
the extension of tidal theory to an Earth-like exoplanet unre- 
liable. Measurements of the martian dissipation function find 
Qmars = 85.58 £0.37 (Bills et al., 2005) or Qmars=79.91 + 0.69 
(Lainey ef al., 2007). The value for Mercury is Qmercury < 190, 
and for Venus it is Qyenus<17 (Goldreich and Soter, 1966). 
Observations of Io imply Q;.<100 (Peale et al., 1979; Segatz 
et al., 1988), while the Moon lies at Qyoon= 26.541 (Dickey 
et al., 1994). These measurements indicate a similar order of 
magnitude for the tidal dissipation function of all desiccated 
bodies in the Solar System. 

Tidal effects on Venus may be particularly complicated due 
to its large atmosphere. The rotation rate of Venus is such that 
it always presents the same face to Earth at inferior conjunc- 
tion. Although this may be by chance due to Venus ap- 
proaching tidal lock with the Sun, it has also been suggested 
that at some point in the past thermal atmospheric tides could 
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have provided an equal but opposite torque than that from the 
rocky component. If that occurred, the torque from Earth could 
dominate and spin-lock Venus to the Earth (Gold and Soter, 
1969). More recent and detailed work has shown that this in- 
teraction can be significant and permits chaotic spin evolution 
that eventually settles into one of four possible states (Correia 
and Laskar, 2001). Although a large atmosphere can play a 
significant role in the evolution of terrestrial planets (Cor- 
reia and Laskar, 2003), we ignore their effect here due to the 
large number of unknowns for an exoplanet. In real bodies, 
Q is a function of the bodies’ rigidity mu, viscosity y, and 
temperature T (Segatz et al., 1988; Fischer and Spohn, 1990). 
A comprehensive tidal model would couple the orbital and 
structural evolution of the bodies since small perturbations 
in T can result in large variations in Q (Segatz et al., 1988; 
Mardling and Lin, 2002; Efroimsky and Lainey, 2007). 
Moreover, tidally generated heat will not necessarily be 
transported to the surface quickly. The internal structure of 
the planet may trap the energy for long periods of time, 
resulting in catastrophic tectonics when the heat is finally 
released, as may have occurred on Venus (Turcotte, 1996). 
This episodic volcanism could change tye, dramatically, but 
on the other hand, complete lithospheric overturn is prob- 
ably just as detrimental to life as desiccation. Nonetheless, 
future work should explore these details. Henning et al. 
(2009) performed detailed calculations and argued that 
ky=0.3 and Q=50 are reasonable choices for dry planets. 

The tidal properties of stars are about as poorly con- 
strained as they are for planets. The primary constraint 
comes from the distribution of orbits of binary stars (e.g., 
Zahn, 1977), whose orbits are circularized by the tidal in- 
teraction. However, for the closest binaries, where the tidal 
effects are strongest and one might think the constraints the 
tightest, the situation is complicated by the early phases of 
stellar evolution (Zahn and Bouchet, 1989), in which stellar 
radii contract rapidly. As the radius enters the tidal evolution 
at the 5" power, the radial contraction is a major effect. 
Khaliullin and Khaliullina (2011) revisited this point and 
found that, although radial contraction does play a major 
role, they were unable to combine models of the equilibrium 
tide and radial contraction to match observations. Weinberg 
et al. (2011) tested the linear theory against nonlinear models 
and found that in general the linear theory does not capture 
the physics of the more sophisticated approach. These 
studies demonstrate that there still exist serious problems 
with the linear theory, and quantitative discrepancies exist. 
However, it does qualitatively model the relevant physics 
and allows wide sweeps of parameter space to be a com- 
putationally tractable problem. We therefore continue to use 
it but explicitly acknowledge that known weaknesses persist. 

A wide range of values for the tidal quality factor and tidal 
time lag for stars is present in the literature, ranging from 
QO,=10* (Adams et al., 2011) all the way to Q,. =10? (Matsu- 
mura et al., 2010), with many others finding intermediate 
values (e.g., Lin et al., 1996; Carone and Patzold, 2007; Jack- 
son et al., 2008b, 2009; Ibgui and Burrows, 2009). For the star, 
we take an intermediate and typical value of Q,=10°. 
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